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ABSTRACT 



We report on very long baseline interferometry (VLBI) observations of the radio emission from the inner region of the Crab Nebula, 
made at 1 .6 GHz and 5 GHz after a recent high-energy flare in this object. The 5 GHz data have provided only upper limits of 0.4 milli- 
Jansky (mjy) on the flux density of the pulsar and 0.4mJy/beam on the brightness of the putative flaring region. The 1.6 GHz data 
have enabled imaging the inner regions of the nebula on scales of up to w 40" . The emission from the inner "wisps" is detected for the 
first time with VLBI observations. A likely radio counterpart (designated "CI") of the putative flaring region observed with Chandra 
and HST is detected in the radio image, with an estimated flux density of 0.5 ± 0.3 mJy and a size of 0''2-0"6. Another compact feature 
("C2") is also detected in the VLBI image closer to the pulsar, with an estimated flux density of 0.4 ± 0.2 mJy and a size smaller than 
0"2. Combined with the broad-band SED of the flare, the radio properties of CI yield a lower limit of « 0.5 mG for the magnetic field 
and a total minimum energy of 1.2 x 10 41 ergs vested in the flare (corresponding to using about 0.2% of the pulsar spin-down power). 
The 1.6 GHz observations provide upper limits for the brightness (0.2 mjy/beam) and total flux density (0.4 mJy) of the optical Knot 
1 located at 0''6 from the pulsar. The absolute position of the Crab pulsar is determined, and an estimate of the pulsar proper motion 
(pi a = -13.0 ± 0.2mas/yr, fi s = +2.9 ±0.1 mas/yr) is obtained. 

Key words. ISM: supernova remnants — radio continuum: stars — pulsars: general — pulsars: individual: Crab pulsar 



1. Introduction 

The inner region of the Crab Nebula contains a complex and 
dynamic structure with a number of elliptical ripples ("wisps") 
varying in flux on ti mescales of days and moving outwar d with 
speeds of up to 0.7 c (iTanvir et al.ll997l : lHester et al.l2002h in the 
optical, while radio meas urements with the VLA yield somewhat 
lower speeds of * 0.3 c dBietenholz et aljr2004l) . The origin of 
the wisps is not yet known, but they are generally thought to be 
associated with instabilities downstrea m of the shock in the win d 
from the pulsar that powers the nebula dGallant & Aronsiri994l) . 

The wind derives its power from the pulsar spin-down inject- 
ing an energy o f ~ 5 x 10 38 erg/s into relativistic particles an d 
magnetic field (iRees & Gunnlll974t iKennel & Coronitilll984l) . 
The particles are thought to be mostly electrons and positrons, 
accelerated to relativistic energies , with possibly a lesser num- 
! ber of ions dArons & Tavanill 1994b . 

The exact mechanism by which the spin-down energy 
is released, transpo rted, and dissip at ed remains poor! 
understood (ct, Ichedia et al.l Il997t iHester et all 1998I 
' Lvutikovl I20031: ISpitkovskv & Aronsl 1200 
with shocks, plasma insta- 
bilities, pair acceleration, Poynting flux, and magnetic field all 
suggested to affect the evolution of the wind. 

Fast variability possibly related to abrupt energy releases 
in the nebula was first reported in the radio on time scales 
of ~ 1 day, corresponding to an emitting region of < O'.'l in 
size dMatvevenkdll975blMatvevenko & Kostenkoiri979l) . On 22 
September 2010, a gamma-ray flare was discovered for the first 
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time from the directio n of the Crab Nebu la with the AGILE pair- 
production telescope dTavani et al.ll201 II) . offering a rare oppor- 
tunity to follow in detail the evolution of a flaring r egion in the 
nebul a. The flare was confirmed with Fermi/LAT dAbdo et al.l 
1201 1[) indicating a variability time-scale as short as 12 hours 
dBalbo et al.ll2011h . 

The flare was not detec ted in the keV regime with Swift/XRT 
dEvangelista et alJ 1201 Oh and INTEGRAL dFerrigno et al.1 
l2010bl) . The Fermi/LAT monitoring had further shown an 
abrupt decrease of the gamma-ray flux on 23 September, with 
the emission returning to its pre-flare level. Based on the 
short duration of the high-ener gy flare, it has been suggested 
dKomissarov & Lvutikovl 1201 Oh that the flaring material is 
located in the optica l Knot 1 separated by » 0'.'6 from the pulsar 
dHester et al.l 1 1995b . However, subseq uent high-resolution 
imaging with Chandra (on 28 September. iTennant et al.ll2010h 
revealed three stationary compact (« 1" in size) regions (la- 
belled A,B, and C) of enhanced emission located on the X-ray 
bright ring. Similar structures are are also visi ble in an HST 
image of the Crab Nebula taken on 02 October dCaraveo et al.l 
120101) . These kno ts present themselv es as a viable site for the 
flaring emission dTavani et al.l 1201 ll) . and their compactness 
may explain naturally the non-detections in observations with 
lower spatial resolution where the extended nebula emission 
dominates. The Chandra and HST observations of the knots 
indicate a strong, rapid, and efficient energy release in the 
flaring material and pose the question of the relation between 
the optical, X-ray, and gamma-ray emission. Given the position 
of the knots within the nebula (located not too far off the jet 
axis and slightly closer than the innermost wisp), it is not clear 
whether the flaring region is associated with the jet or with the 
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Table 1. Antenna parameters 



Antenna 


D 


^obs 


SEFD 


"min "max 




[m] 


[cm] 


[Jy] 


[km] 


Effelsberg 


100 


6,18 


20,19 


266- 853 


Westerbork 


66° 


6,18 


60,30 


266-1002 


Jodrell Bank 


25 


6,18 


320,320 


17-1388 


Medicina 


32 


6,18 


170,600 


757-1429 


Onsala 


25 


6,18 


600,320 


601-1429 


Torun 


32 


6,18 


220,230 


637-1441 


Yebes 


40 


6 


90 


1313-2152 


Cambridge* 


32 


6 


136 


197-1313 


Darnhall* 


25 


18 


356 


17-1404 


Knockin 6 


25 


18 


356 


67-1441 


Hartebeesthoeck' 


26 


18 


450 


7453-8525 



Notes: a - for Westerbork operating as a phased array; b - MERLIN 
antennas; c - Hartebeesthoeck was added ad hoc in the observation 
at 18 cm, while participating in the e-VLBI session on November 5. 
Column designation: D - antenna diameter; A obs - observing wave- 
length used; SEFD - system effective flux density describing the an- 
tenna sensitivity at corresponding wavelength; S m i„-fi max - range of the 
telescope separations from a given antenna (for European telescopes, 
baselines to the South African antenna in Hartebeesthoeck are not con- 
sidered). 

Table 2. Target and calibrator sources 



Object/Coordinates 



o- pos [mas] 





Crab-A" 




CJ2000 


05 h 34 m 32?332 


-10 


<5|2000 


+22° 00' 5Z'45 
J0518+2054 6 
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fnooo 
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Fig. 1. Spatial frequency coverage (Fourier domain or uv- 
coverage) of the VLBI visibility data for the target source at 
1 .6 GHz. The coverage is plotted in Megawavelength (MA) units, 
referring to spatial frequencies of fi/ jm /(/t bs X 10 6 ), where B[ m 
is the projected baseline length for a given pair of antennas 
(/, m). Baselines to Hartebeesthoeck, covering a range of 36- 
48.8 MA, are not shown. The smallest spatial frequency sam- 
pled is 10~ 2 M/i, corresponding to »2" , but the central gap in 
the Mv-coverage results in the structural sensitivity substantially 
reduced for scales larger than 0"2. 



Notes: <x pos - position errors in right ascension and declination. 
References: a - iTavani et al.l201ll and this work; b - iPetrov et alj 2008. 



equatorial outflow. Deciding on these issues relies on detailed 
imaging of the emitting region, which is best achieved in the 
radio regime. 

We present here high-resolution wide-field radio images of 
the flaring knot and the central region of the nebula, obtained 
with very long baseline interferometry (VLBI) observations at a 
wavelength of 18 cm. These observations probe a wide range of 
angular scales, from ~ 10 milliarcseconds to ~ 40 arcseconds, 
detect the flaring region, and reveal, for the first time, the intri- 
cate morphological structure of the radio emission within about 
20 parsecs of the pulsar. The observation and data reduction are 
described in Section 2. The resulting images are presented in 
Section 3, and discussed in the context of physics and evolution 
of the flare in the Crab Nebula. 



2. Observations and data analysis 

The flaring region was observed with VLBI on 5 November 20 10 
at 1.6GHz(/l obs = 18 cm) and on 23 November 2010 at 4.9 GHz 
(/lobs = 6 cm). The total of eight EVl\j telescopes and three 



MERLINr[telescopes participated in the observations, operating 
in e-VLBJj mode. Basic parameters of the participating anten- 
nas and ranges of the respective baseline lengths, B, are given in 
TableQ] 

The 1 .6 GHz observations were made over 1 1 hours at a cen- 
tral reference frequency of v b s = 1616.49MHz. The 4.9GHz 
observations were made over 8.2 hours at a central reference 
frequency of v ODS = 4948.74MHz. The 6 cm and 18 cm obser- 
vations have sampled ranges of visibility (uv) spacings, B/A, of 
2-36 x 10 6 A (2-36 M^) and 0.01-49 MA, respectively (note that 
Hartebeesthoeck took part only in the 18 cm observation). 

For both observations, the field of vie w was centered o n 
the position of the flaring region, Crab-A dTavani et alj|201 11) . 
as measured from the HST data. The nearby bright (S i.6GHz = 
0.42 Jy) and compact radio source J05 1 8+2054, separated by 4!0 
from Crab-A, was used as a phase-reference calibrator. A cali- 
brator/target cycle of 2/5 minutes was used, with about 30 sec- 
onds allocated for slewing. Two other calibrators, J0530+1331 
and 3C48, were observed as fringe finders. Positions and flux 
densities of the target and phase-reference calibrator are given 
in Table [2] 



European VLBI Network, www.evlbi.org 



- Multi Element Radio Linked Interferometer Network, 
www.merlin.ac.uk 

3 electronic- VLBI, a mode of VLBI operations, with direct links be- 
tween participating telescopes and the correlator maintained via optical 
fibre connections; http://services.jive.nl/evlbi/ 
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The EVN data were recorded at a rate of 1024Megabit/s 
(Mbps), in dual-circular polarisation, which yielded and observ- 
ing bandwidth of 1 28 MHz per polarisation. The total bandwidth 
was divided into 8 intermediate frequency (IF) bands, each cov- 
ering 16 MHz and split in 32 spectral channels of 0.5 MHz in 
width. The MERLIN data were obtained at a 128 Mbps recoding 
rate, in dual-circular polarisation, yielding an observing band- 
width of 16 MHz per polarisation. The MERLIN data, corre- 
sponding to a single IF band of the EVN data, were recorded 
also using 32 spectral channels of 0.5 MHz in width, for each po- 
larisation channel. The recording was made in a two-bit record- 
ing mode at all participating antennas. The resulting coverage of 
the Fourier domain (wv-coverage) provided by the visibility data 
on the target source is shown in Fig. [1] for the observations at 
1.6 GHz. 

Correlation of the data was done at the EVN correlator facil- 
ity of the Joint Institute for VLBI in Europe (JIVE). The output 
correlated data were produced with a 2 seconds averaging time. 
The combined choice of the integration time and the frequency 
channel width ensured that a field of view as large as 1'6 (at 
1.6 GHz) and 0'5 (at 5 GHz) could be imaged using only the 
visibility data from the European baselines (B < 10 MA) shown 
in Fig.Q] with the synthesised beam degraded by less than 10%. 
No pulsar gating was applied (this correlation mode is currently 
unavailable for e-VLBI observations). 

The correlated data were processed using AIP5J3 software. 
The data scans were first flagged using automatic flag tables pro- 
duced during the correlation. The visibility amplitudes were cal- 
ibrated, based on the system temperature measurements and sta- 
tion gain information. The data for the phase-reference calibrator 
were then fringe-fitted, and a bandpass calibration was applied. 
The phase-reference calibrator was imaged and found to have 
only a weak extended structure, ensuring that it was not con- 
taminating the phase solutions. The phase solutions were then 
transferred to the target scans, and phase-calibrated data for the 
target source were imaged. 

2.1. 5 GHz data 

No positive detection of emission can be made from the 5 GHz 
data for the flaring region and for the pulsar itself, even in a 
heavily tapered image with a beam (FWHM) of 0'.'12 x 0"03. 
The pulsar emission is reported to have a continuum flux den- 
sity of 14.4 ± 3.2 mJy at 1.4 GHz and a very soft spectral index 
o-ps, = -3.1 + 0.2 (lLorimer et al.lll995l) . This corresponds to a 
flux density of a 0.3 mJy at 5 GHz. We estimate an r.m.s. noise 
of cr 6cm ss 0.16mJy/beam, from the inverted 5 GHz visibility 
data. Thus the pulsar flux is expected to reach an « 2 <x level in 
the 5 GHz image. Based on these considerations, an upper limit 
of ~ 0.4 mJy can be provided for the flux density of the pulsar. 
The corresponding upper limit on the brightness of the knot A is 
then ~ 0.4 mJy/beam, for the above mentioned beam size. 

The following discussion is therefore focused primarily on 
the results obtained from the VLBI data at 1.6 GHz. 



2.2. Effects of the pulsar and extended nebula on imaging 

The presence of the pulsar and an extremely large (420" x290" ) 
and bright nebula results in an increase in the image noise level 
in the flaring region. For the antenna configuration used in the 
observation, sensitivity reductions by factors of 3.0 and 3.2 are 



4 Astronomical Image Processing Software, National Radio 
Astronomy Observatory, USA 



expected for the EVN and MERLIN parts of the data, respec- 
tively. These estimates account for the 1.6 GHz flux density 
(~ 900 Jy) of the nebula and the primary beams response of the 
EVN and MERLIN antennas to the nebula. The resulting ex- 
pected r.m.s. image noise in the combined data is as 60juJy/beam 
(in a full-resolution image made from naturally-weighted visi- 
bility data). The contribution of the pulsar to the image noise is 
not as significant. At the observing frequency of 1.6 GHz, the 
expected continuum flux of 8.6 mJy, and sidelobes from the pul- 
sar resulting from deconvolution defects are expected to become 
relevant at a ~ 10 /iJy/beam brightness level, which is lower than 
the estimated thermal noise. 



2.3. Resolution and structural sensitivity 

The visibility data at 18 cm sample a range of t/v-spacings from 
10~ 2 MA to 48.8 MA, with a large gap between 9 MA and 36 MA. 
The full-resolution, image made from the uniformly-weighted 
data, including the baselines to Hartebeesthoeck, has a synthe- 
sised beam with a full width at half maximum (FWHM) of 
19.9 x 2.2 mas, while the European baselines yield a beam with 
a FWHM of 25.7 x 22. 1 mas. The combination of angular scales 
sampled and flux density sensitivity of the data on the baselines 
to Hartebeesthoeck makes them useful only for imaging and as- 
trometry of the pulsar itself. 

The largest formally detectable scale in the data is about 20" 
(«20% of the unaberrated field of view). However, this scale is 
sampled only by a single baseline, and a realistic limit on ad- 
equate structural sampling is a;0"2 («1 MA), noting the gap in 
the Mv-coverage shown in Fig. Q] This gap corresponds to an- 
gular scales of 0"2-2" . It should be thus expected that smooth 
emission on scales larger than 0"2 (and in the 0"2-2" range in 
particular) could only be poorly recovered from the data (lead- 
ing to "spotty" structures in images obtained using deconvolu- 
tion, or requiring progressively higher signal-to-noise ratio on 
these scales in order to be successfully recovered by the maxi- 
mum entropy method). The limited structural sensitivity of the 
VLBI data acts effectively as a low frequency uv-filter, reducing 
the contribution from large spatial scales where diffuse emission 
from the extended nebula dominates strongly. This leads to a bet- 
ter detectability of bright and compact regions which otherwise 
would have been swamped by the bright diffuse emission. 



2.4. Imaging 

Taking into account the issues of resolution and structural sensi- 
tivity outlined above, we have made three different images from 
the data: 1) a full-resolution image of a small region around the 
pulsar (Fig. |2); 2) a single-scale CLEAN image aiming at de- 
tecting the flaring region (Fig. [3); and 3) a multi-scale CLEAN 
image of extended emission within a region of about 40" in size 
(Fig. |2). The latter two images have been made using only the 
data from European baselines covering the wv-spacings of up to 
10M/1 (this has been done in order to improve the sensitivity 
to extended, low surface brightness emission). Properties of the 
images obtained are summarised in Table [3] 



2.4.1. Pulsar 

The full-resolution image of the Crab pulsar is made from the 
entire data, using the uniform weighting and restricting the field 
of view to 0"3. A phase shift of Aor = -5"479 in right ascension 
and A<5 = -0"273 is applied to the data in order to center the 
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Table 3. Properties of eVLBI images of the Crab Nebula at 
1.6 GHz 



Image 


Beam 
[arcsec] 


^ peak 

[mjy/beam] 


0"rms 

[mJy/beam] 


1 (Fig. 2) 

2 (Fig. 3) 

3 (Fig. 4) 


0.020 x 0.002 
0.150x0.150 
0.500 x 0.500 


8.35 
8.18 
8.53 


0.28 
0.20 
0.16 




2 


4 


6 



22 00 52.24 



52.22 



O 

l- 
< 




52.12 



05 34 31.942 31.941 31.940 



31.939 31.938 31.937 31.936 31.935 31.934 
RIGHT ASCENSION (J2000) 



Fig. 2. Full resolution image of the Crab pulsar obtained from 
uniformly weighted data. The image has a peak flux density of 
8.3mJy/beam and an rms noise of 0.3mJy/beam. The synthe- 
sised (restoring) beam is 19.9x2.2 mas, oriented at a P.A. = 83° . 
The extreme ellipticity of the restoring beam results from a very 
limited East- West Mv-coverage on baselines > lOM/t. 



image on the location of the pulsar. The resulting high-resolution 
phase-referenced image of the pulsar is shown in Fig. [2] 

Properties of the pulsar emission are estimated by fitting 
an elliptical Gaussian to the image brightness distribution. This 
yields an integrated flux density of 8.4 + 0.5 mJy, which is 
cl ose to the valu e estimated from the flux density measurements 
of lLorimer et alj (1995). Deconvolution of the fitted extent of 
the Gaussian indicates that the pulsar emission is not resolved, 
with the limits on its size obtained from the uncertainties of the 
fit cr ma j = 0.40 mas (at a position angle of 81°) and cr m i n - 
0.09 mas for the major and minor axis of the Gaussian, respec- 
tively. A similar estimate (cr ma j = 0.40 mas, cr m \ n - 0.04 mas) 
is obtained using the SNR-based resolution criterion dLobanovl 
12001 . 



2.4.2. Image of the flaring region 

As described above, extended emission of the Crab nebula has 
been imaged with the data from the European baselines (B < 
10 MA), which has been done in order to improve the sensitivity 
to emission on scales > 20 mas resolved out on the baselines to 
Hartebeesthoeck. 



52 
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Fig. 3. Image of the pulsar and the flaring region. The image is 
obtained by applying a single-scale CLEAN deconvolution on 
the Mv-tapered visibility data on the European baselines. The pa- 
rameters of the image are given in Table [3] For the presentation 
purpose the image is convolved with a 0'.T5 beam (three times 
the natural restoring beam). 



The visibility data are phase-shifted in right ascension by 
Aq, = -6"5, in order to fit better the field of interest within a rect- 
angular grid. The grid extends by 20"9 from the phase-center, 
thus covering almost the entire unaberrated field of view (allow- 
ing for about 5% reduction of the peak response due to band- 
width and time average smearing). As the flaring region is likely 
to be extended, the visibility data are weighted using the natu- 
ral weighting and tapered using a Gaussian taper with <x taper = 
2.2 MA. This yields a synthesised beam of 54 mas x 49 mas, 
which is about 2.5 times larger than the beam obtained without 
the taper. 

The image of the flaring region shown in Fig. [3] is obtained 
by using single-scale CLEAN deconvolution and convolving the 
results with a larger beam (0"15) in order to enhance the con- 
trast for the regions of weak, extended emission. The likely radio 
counterpart of the flaring component (denoted as CI) is clearly 
detected in this image, albeit at a somewhat smaller separation 
from the pulsar than the optical/X-ray emitting region Crab-A. 
The weak radio jet (located in the direction of the "sprite" iden- 
tified in the HST images; cf. lHester et al.ll2002t iMelatos et all 
l2005h and another, weaker but more compact component (C2) 
are also visible in the image. The possible nature of the two com- 
pact components and their relation to the flare will be discussed 
below. 

The bright optical knot (Knot 1) seen in the HST images 
at a 0"6 distance from the pulsar (cf., iHester et al.|[l995l) is not 
detected in the 1.6 GHz image. An upper limit on its brightness 
is provided by the rms (0.2mJy/beam) measured in the region 
where it is located. The non-detection of Knot 1 may imply both 
its weakness in the radio and large extension resulting in it being 
resol ved out in the VLBI data. For the Knot 1 dimensions (0"50x 
0'/16. lHester et al.lll995l) . the upper limit on its total flux density 
at 1.6 GHz is * 0.4 mJy. 
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Fig. 4. Image of the central region of the Crab nebula obtained from naturally weighted data, Mv-tapered data on baselines up to 
10M/L The image is obtained using a multi-scale CLEAN deconvolution and restored with a circular beam of 0"5 in size. The 
image has a peak flux density of 8.5 mJy/beam and an rms noise of 0. 16 mJy/beam. The total flux density recovered in the image is 
148 mJy. 



2.5. Fidelity of detection of the compact features 

The compact components CI and C2 are strongly resolved, and 
they cannot be detected without tapering of the visibility data. 
This implies that their sizes should be a priori larger than ~ 
40 mas and that both emitting regions most likely have a smooth 
extended structure that can be adequately imaged only with a 
dense sampling on short baselines (B < 250 km). Because of the 
lack of such baselines in our data, the appearance and relative 
brightness of CI and C2 depend on the choice of deconvolution 
and image restoration procedures. In the single-scale CLEAN 
image in Fig. [3] which is biased to higher spatial frequencies 
(hence smaller angular scales), the peak brightness is higher in 
the more compact feature C2. The multi-scale CLEAN image 
in Fig. [4] accounts better for lower spatial frequencies (larger 



angular scales), yielding a larger flux density for the more ex- 
tended feature CI, with some of this emission possibly coming 
from the extended emission of the nearby "wisp". Considering 
all these factors, and taking into account the likely angular extent 
of both components, we estimate that the application of single- 
scale CLEAN deconvolution yields detections at an SNR of 4.2 
and 3.2, for CI and C2 respectively. The respective SNR values 
are reduced to 3.8 and 2.5 in the multi-scale CLEAN image in 
which the cumulative spatial sensitivity is shifted to scales that 
are much larger than the angular extent of CI and C2. 

To test the fidelity of detection of such features in the EVN 
data, we amend the visibility data with two sets of simulated 
Gaussian components representing CI and C2 at different posi- 
tion angles and distances from the pulsar. The representations of 
CI have a flux density of 0.5 mJy and a size of 0"2; the respective 
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values for C2 are 0.4 mJy and 0"15. Random noise is added to 
the visibility data, to accommodate for the total signal increase 
due to the simulated features. As the sum flux density of the sim- 
ulated components (1.8 mJy) represents only a small fraction of 
the total flux density in the image (140mJy), the resulting noise 
factor is also relatively small and leads to an increase of the av- 
erage visibility noise from 73.2 mJy to 74.1 mJy. 

Locations of the first set (CI a, C2a) of the added components 
are obtained by rotation of the positions of CI and C2 by -122° . 
For the second set (Clb, C2b), a rotation by -47° is combined 
with a factor of 4 stretch in distance from the pulsar. This en- 
sures that the resulting locations of the simulated components 
span a range of pulsar separations, Ar, ranging from 1"6 to 16"4 
and fall outside of the areas with detected emission. The mod- 
ified data are then imaged using the same procedure as applied 
for obtaining the single-scale CLEAN image. All four simulated 
features are detected in the resulting image and the ratios, (s, of 
their fluxes to the fluxes of the original components are close to 
unity (see Table [H, and only the most distant simulated feature 
shows a perceivable decrease in flux density, owing most likely 
to the smearing effects. These results support further the reliabil- 
ity of the detection of the compact features CI and C2. 



Table 4. Flux density ratios for simulated compact features 



Name 


Ar 


P.A. 


& 


C2a 


L'6 


-49? 


1.11 ±0.51 


Cla 


A'!\ 


-41°0 


0.91 ±0.32 


C2b 


6'.'6 


26!0 


0.96 ± 0.47 


Cla 


16:4 


34!0 


0.79 ±0.31 



Notes: Ar, P.A. - component distance and position angle with respect 
to the pulsar; f s - measured ratio between the flux density obtained in 
the test image for the simulated and the original feature. 



2.5.1 . Image of the central region of the nebula 

Beyond the extent of the flaring region, the CLEAN algorithm 
finds positive flux as well, indicating likely detection of the 
two inner "wisps" visible in the optical and rad io images ob- 
tained with HST and VLA dBietenholz et al.l2004 . As expected, 
this emission is not well recovered by CLEAN, resulting in a 
"spotty" appearance of the elliptically-shaped wisps. In order to 
improve the representation of large-scale structures in the image, 
we have applied multi-scale CLEAN (MS-CLEAN) algorithm. 

Application of MS-CLEAN leads to moderate improvement 
of the image at large scales (Fig. @}, and it has also resulted in 
a moderate reduction of the r.m.s. noise in the image. The dis- 
tribution of the residual flux in the image is shown in Fig. 
The distribution is essentially Gaussian, implying an r.m.s. of 
0.16mJy/beam. At the same time, the application of multi-scale 
CLEAN has introduced a substantial power on scales > 0"5 and 
thus reducing the contrast on smaller scales (which is visible in 
particularly in Fig. |4]at the location of the weaker feature C2). 

The total CLEAN flux density in the image is 148 mJy, 
which reflects only a fraction of the total flux density in this 
area (as the VLBI data do not sample spatial scales larger than 
20" and indeed have problems with recovering flux on spatial 
scales larger than about 0"2. 
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Fig. 5. Histogram of the residual flux distribution in the VLBI 
image, after applying the multi-scale CLEAN. The histogram is 
well approximated by the Gaussian noise (solid line) with an 
r.m.s. of 0.16mJy/beam. 



2.5.2. Fidelity of the large-scale structure 

The large-scale structure detected in the image shown in 
Fig. Bean be compa red to the structure imaged with the VLA 
(Bietenho lz et al.l | 20 04). This comparison is presented in Fig [6] 
which relates the VLA image to the multi-scale CLEAN image 
restored with a 1'.'4 beam (to match the VLA image resolution). 
The location and shape of the inner "wisps" agree well in the 
VLBI and VLA images. It should be noted of course that the 
two images are not contemporary, and thus the apparent agree- 
ment may be simply fortuitous. Nevertheless, this argues for the 
overall plausibility of the larger scale structure in the VLBI im- 
age. 

One particularly remarkable difference between the structure 
seen in the VLA and VLBI images is observed in the south- 
western direction, where the VLBI image does not show any en- 
hancement in flux, while a strong and extended emission region 
is visible in the VLA image. It is presently impossible to judge 
whether this disagreement is caused by the evolution of the ra- 
dio emission or by deficiencies in the VLBI data and imaging. It 
should be noted however that this region is situated at more than 
30" (^600beamwidths) away from the phase-center of the VLBI 
data, and the resulting reduction in peak response and beam de- 
terioration may reduce detectability of weak and extended emis- 
sion. 

The imaging results summarised above indicate that the 
VLBI observations have not only been sensitive enough to de- 
tect the presence of a compact flaring emission in the Crab 
Nebula, but they also enable obtaining a detectable response 
from the inner wisps of the nebula. However, the fidelity of lo- 
calisation of this emission remains marginal, given the limited 
wv-coverage of the VLBI data. It is obvious that improved uv- 
coverage on shorter baselines is required in order to be able to 
image these structures, and a full-fledged combination of eMER- 
LIN and EVN, both operating at 1 Gbit/s recording rates, would 
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Fig. 6. The VLBI image of the Crab nebula (grey scale) overlaid 
with a VLA im age (contours) of the same region obtained by 
iBietenholz et al . (2004). The VLBI image is obtained from the 
multi-scale CLEAN image by removing the contribution from 
the pulsar and restoring the resulting image with a 1'.'4 beam in 
order to match the resolution of the VLA image. Locations of 
the inner "wisps" agree well in the two images, despite potential 
complications due to the non-contemporary observing epochs 
and possible image artifacts in the VLBI image. See further dis- 
cussion in the text. 



provide substantially better images of the central region in the 
Crab nebula. 

The measured brightness of the wisps at 1 .6 GHz and the 
r.m.s. estimate at 5 GHz can be reconciled with the "canoni- 
cal" spectral index of a = -0.3 obtained i n the radio r egime 
for the total flux de nsity of the nebula dBaars et al.l Il977t 
iKovalenko et al.lfl994l) as well as for the reso lved structures in- 
cluding the inner wisps ( IBietenholz et al.l 19971) . The inner wisps 
detected in the VLBI image at 1.6 GHz have an average bright- 
ness of 0.3 mJy/beam. Adopting a = -0.3 yields an expected 
5 GHz brightness of 0.2 mJy/beam, which compares well with 
the r.m.s. of 0. 16 mJy/beam estimated from the 5 GHz data. 



ered in the calibrator image is 290. 1 + 1.8 mJy, and the peak flux 
density is 271.0 + l.OmJy/beam, further underlying the com- 
pactness of the structure. The calibrator image has a convolving 
beam of 18 x 12 mas and an r.m.s. noise of 0.8mJy/beam (and 
the respective peak-to-noise SNR of 340), thus limiting the an- 
gular size of the calibrator to < 1.1 mas and introducing an error 
of 0.06 mas to the positional measurements (lLobanovl2005l) . For 
the calibrator, we obtain the following J2000 position (the errors 
listed are formal errors of the fit): 

ajiooo = 05:18:03.8245014+0.0000014, 



3J2000 



= 20:54:52.497662 + 0.000028 . 



Analysis of the high-resolution image of the pulsar yields an ap- 
parent J2000 position (the errors are formal errors of the fit): 

ajaooo = 05:34:31.9383014 + 0.0000081 , 

dnooo = 22:00:52.17577 + 0.00018 . 

Taking into account the positional errors of the calibrator 
listed in Table [2] more conservative estimates of the errors of 
the pulsar position are given by A a = 0.18 mas (0?000012) 
and Ag = 0.22 mas. These estimates do not include errors 
due tropospheric and ionospheric delays and errors from the 
calibrator-target phase transfer normally considered in measure- 
ments made with dedicated VLBI astromet ry observations (cf., 
iGuirado et all 1 9951: iMartt-Vidal et al.ll2~0 08). Based on the anal- 



ysis of these factors (IGuirado et al.l 1 19951) . we can expect that 



their contribution to the positional errors should be within « 
0.05 FWHM [mas] A0[deg] milliarcseconds, with A9 giving the 
calibrator-target distance. For the FWHM of the full resolution 
image, one can therefore expect that the positional errors can 
reach Aa = 3.4 mas and A6 - 0.45 mas (reflecting a much better 
north-south resolution of the image). 

The position measurement made from the VLBI data can be 
combin ed with the Crab pulsar positions obt ained from pulsar 
timing dTavlor et al.lll993J lHan & Tianll 19991) . yielding the fol- 
lowing estimate of the pulsar proper motion: 

fi a = -13.0 + 0.2 mas/yr, 

Us - +2.9 + 0.1 mas/yr. 

This estimate agrees well wit h recent estimates made from a long 
range of HST observations dNg & Roman j|2006t iKaplan et al.l 

2008). 



3. Results and Discussion 

The VLBI observations of the Crab nebula have yielded accurate 
positional and photometric information about the pulsar itself 
and the flaring region, as well as marginal morphological and 
photometric information about the inner wisps of the nebula. 



3.1. Astrometry and proper motion of the pulsar 

The position of the pulsar is measured relative to the position 
of the phase-reference calibrator, J0518+2054. The position of 
J0518+2054 has been determined with an accuracy of 0. 12 mas, 
both in right ascension and declination dPetrov et al.l l2008). The 
calibrator image obtained from the VLBI data shows no signif- 
icant structure that could affect the calibrator phases and intro- 
duce additional positional errors. The total flux density recov- 



3.2. Flaring region 

The flaring components CI and C2 are both very weak, detected 
at SNR < 4, which makes it difficult to provide robust estimates 
of their properties. Table [5] lists our best estimates of the posi- 
tions, flux densities, and sizes of the two components, based on 
Gaussian fits to the single-scale and multi-scale image (with the 
latter providing a somewhat better account of the extent of the 
emitting regions). We use these positions to estimate brightness 
temperature of the emission in respective regions. The resulting 
values indicate brightness temperature in excess of « 4000 K (it 
should be noted here that estimating such low values of bright- 
ness temperature is enabled by the presence of short baselines 
provided by the MERLIN antennas). The upper limit on the size 
of CI corresponds to the maximum size of a 0.5 mJy feature that 
could be detected in the data. 
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Fig. 7. Image of the radio emission (green contours; multi-scale CLEAN image, convolved with a 0"7 be am; contours start a t 
SNR — 2 and increase with an interval of 0.5) overlaid with ima ges of the optical (left panel, 02 October 2010: ICaraveo et aDl2010h 
and X-ray (right panel, 28 October 2010; Ferrigno et al.ll2010ah emission in the flaring region. The radio emission from the flaring 
region CI is offset from the knot A dTavani et alj|201 lb identified in the optical and X-ray emission. The overall appearance of the 
optical and radio emission is similar in the flaring, with a clear "knee-like" structure seen in the optical and trailed by a similar 
structure in the radio. It should also be noted that the apparent discrepancy in appearance of CI and C2 in this figure and in the 
image in Fig. 3 (convolved with a 0"15 beam) results from the fact that the component C2 is compact (< 0"2 in size) while the size 
of CI can be as large as 0"6 (see discussion in Sect. 13. 2t 



Table 5. Properties of compact flaring components 





Cl 


C2 


R.A. [s] 


32.2271 ±0.0011 


32.03851 ±0.00052 


Dec. [" ] 


53.096 ± 0.014 


52.8600 ± 0.0070 


S i.6 GHz [mJy] 


0.5 ±0.3 


0.4 ± 0.2 


ei"] 


o:'2 - o% 


<o:'2 


r b [K] 


1000 - 5000 


>4000 



Notes: Si.ghz - component flux density; 6 - estimated size (C2 is un- 
resolved, and a limit on the size is estimated from the SNR following 
lLobanovl II2005J): for Cl, the upper limit on the size is given by the de- 
tection limit in the data); T b - brightness temperature calculated for the 
respective size and flux estimates (note that short baselines of MERLIN 
(10~ 2 -10 _1 M/1) provide sensitivity to emission with brightness temper- 
ature greater than » 1000 K). 



An overlay of the radio imag e of the flaring regio n with 
an HST image from 2 Octobe r dCaraveo et al] 1201 Oh and a 
Chandra image from 28 October dFerrigno et al.l2010al) is shown 
in Fig. [7] The radio image is convolved with a 0'.'7 beam to fur- 
ther emphasise weak extended emission. This overlay shows that 
the radio emission traces well the inner wisp and "sprite" (region 
B), commonly identified in the optical/X-ray images. The radio 
component Cl is located close to a "knee-like" region (region 
A) detected in the optical and X-ray images and believed to be 
associated with the flare dTavani et al.ll201 II) . This morphology 
may also be reflected by the radio emission, but establishing this 
relation clearly requires better quality of radio imaging on arc- 
second scales. 



3.3. Origin of the radio knots 

A significant offset of both Cl and C2 from the jet axis is ev- 
ident. If these two features are associated with the jet (and are 
not part of the wisps produced in the equatorial outflow), it may 
be speculated that they are located on a "wall" formed by inter- 
action of the rotating outflow with the ambient material in the 
nebula. 

The offsets Ar of Cl and C2 with respect to the jet axis, 
Z, directed at a position angle if/j et as 130° (estimated from 
the image in Fig. [3]l can be represented by a power-law ex- 
pansion Arjx^g^^which yields m « 0.7. A jet viewing angle 
of 60° dCheng et al.ll2000t iKomissarov & Lvutikovll2010h is as- 
sumed in these calculations. The expansion derived is consis- 
tent with the one ex pected for the jet boundary in a ma gneti- 
cally confined flow (JFendtl 1 1 9971 iFendt & M emola 2001) with 
a weakly evolving poloidal magnetic field component. One can 
assume that the jet boundary is de fined by the outermost flux 
surface (cf.. IFendt & Mem ola 2001) and that the collimation of 
the flow begins near the light cylinder, R\ c , with an initial radius 
of the flow of w 10^i c . It would then require an electric current 
I p « 4 x 10 14 A, a magnetic field B ^ 4 x I0 n G, and a magnetic 
flux M^p « 6 x 10 24 G cm 2 (all agreeing well with e arlier esti mates 
made for the Crab pulsar; cf.. lHan kins & Eilek 2007; Hirotani 
l2006h to explain the locations of the knots Cl and C2, assuming 
that these knots are located near the jet boundary. We note that 
the value of I p is larger than the Goldreich- Julian current, which 
can be explained by pair cascades yielding particle densities wel l 
in excess of the Goldreich- Julian density (Arendt & Eilek 2002). 
With the jet expansion derived, the jet should have local half- 
opening angles of cpci = 27° and 0c2 = 37° at the locations of 
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Fi g. 8. The spectral e nergy distribution of the nebula and the flaring component CI. The data on the nebula have been compiled 



Mever et all d2010h . Details on the modelling of the synchrotron nebula are provided in that reference as well. In addition to the 



nebula, the SED of the feature CI is shown together with a time-dependent model of the flare (solid red line), with spectra calculated 
at the epoch of the flare maximum (ff ) and one day before (ff-id) and after (ff+id) the maximum 



the radio knots. For the HST Knot 1, the respective half-opening 
angle is 0hst l — 56° , which is close to the val ue of 60° expected 
for this region dKomissarov & Lvutikovl2010h . These arguments 
lend further support to the suggestion that the radio knots are 
connected to the outflow. However, providing a firm conclusion 
on this matter would certainly require a closer followup of an- 
other prominent flare in the Crab nebula. 

3.4. Connection to the high-energy flare 

Based on the observational results summarised above, the faint 
and compact radio features CI and C2 could in principle be re- 
lated to the gamma-ray flare, even if they are not co-spatial with 
the active regions observed in the optical and X-ray bands. In 
the following, we consider a scenario in which the brighter fea- 
ture CI is related to the gamma-ray flare and discuss the conse- 
quences. 

The spectral energy distribution (SED) of CI is shown in 
Fig. [8] together with the nebula emission (see lMever et al.ll2010l 
for details of the data). This SED combines the flux density mea- 
sured at 1.6 GHz, the upper limit obtained at 5 GHz, and the X- 
ray flux of the knot A treated here as a (soft) upper limit on the 
X-ray emission from the immediate vicinity of CI. Assuming 
that the underlying particle distribution present in the region CI 
has the same spectral slope (dn/dy = nay ) as the electrons 
responsible for radio emission in the nebula, the resulting spec- 
trum matches quite well both the radio as well as the gamma-ray 
emission across 14 orders of magnitude in energy. 

The observed variability time-scale of about 24 hrs 
(Bal bo et aD 1201 II) constrains the size of the emitting region 
to be smaller than r < t v . dr c - 3 x 10 15 cm(f var /10 5 sec) = 
0"1 (fvar/10 5 sec). This is consistent with the observed extension 
of CI. Assuming that the gyro-radius r g < cf var and that gamma- 
rays are produced via synchrotron emission, a lower limit on the 
B-fieldissetbyfi>0.5mG(£/100MeV) 1/3 (f var /100ksec)- 2/3 . 



A conservative upper limit on the magnetic field of » lOOmG 
follows from the argument that the total energy in the mag- 
netic field should not exceed the spin down power of the pul- 
sar integrated over the duration of the flare. The equiparti- 
tion magnetic field (assuming minimum total energy) is B eq = 
3 mG(f var /100ksecr 6/7 (£/100MeV) 1/7 . The corresponding to- 
tal minimum energy (w tot = w e + wb = 1.2 x 10 41 ergs). The 
injection power required to generate the burst isw* w>tot/? va r = 
1.2xl0 36 ergs/s = 0.2 % of the available spin-down power. Even 
though this is a comfortable margin, deviations from the min- 
imum energy requirement are limited to roughly one order of 
magnitude. 

Given the limited knowledge about the dynamics of the fea- 
ture CI, it is difficult to estimate directly any effects related 
to relativistic beaming. A kinematic estimate of the Doppler 
factor D = y (1 - yScosi?) by adopting a jet viewing angle 



§ = 30° favoured in several recent works (e.g.. IN g & Romani 
2004 ; iKomissarov & Lvubarskvll2004t iKomissarov & Lvutikov 



20101) and assuming that the jet material moves at a speed sim- 



ilar to the range of speeds observed in the wisps. This yields 
Doppler factors in a 1.4-1.8 range. An upper limit on the D can 
also be derived based upon the non-observation of an accompa- 
nying flare at very high energy gamma-rays. The con tempora- 
neous observatio ns using the VERITAS and MAGIC (iMariottil 
l2010tlOngl2010l) air Cherenkov telescopes rule out any substan- 
tial deviations of the flux in the energy regime between 100 GeV 
and TeV. Taking the relative flux errors of roughly 10 % as an up- 
per limit on any additional flaring component, we can estimate 
an upper limit on£)< (9(100). The actual upper limit depends 
critically on the geometry. For the sake being conservative, we 
assume a synchrotron self-Compton type scenario neglecting ad- 
ditional photon fields. 

These considerations (spectrum and size) are consistent with 
the hypothesis that the region CI detected in the radio image is 
connected with the high-energy flare in the Crab nebula. The un- 
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usually long duration of higher flux state during this event may 
receive further support from the detection of C2, which could be 
a second plasma condensation ejected from the pulsar after the 
main flare in September 2010 and leading to emission enhance- 
ment in the high-energy regime. 

The observed positional offset between the features CI and A 
and the lack of obvious optical/X-ray counterparts for the feature 
C2 raise general questions of the relation between the emission 
detected in all three bands and the nature of the flaring activity in 
the Crab nebula. In the absence of better quality data, these ques- 
tions remain open. As a matter of speculation, one can suggest 
that the flaring activity may be occurring in a relatively thin layer 
close to the jet boundary, thus highlighting the jet edges where 
the path length through the emitting material is larger (indeed, 
the locations of the features A and C can also be reconciled with 
the same jet expansion as derived above for C 1 and C2). Such an 
activity may originate, for instance, from the magnetic reconnec- 
tion or interaction of the jet with the ambient medium. However, 
in the absence of more detailed observational data, it is clearly 
not feasible to provide a reliable judgement on this matter. 

In this respect, we also would like to emphasise that further 
improvement of the quality of radio imaging on arcsecond scales 
(which should be expected after the e-MERLIN becomes fully 
operational) would provide much better reliability for imaging 
of such extended and complicated objects as the central regions 
of the Crab nebula. Therefore, engaging in a detailed monitoring 
program of the Crab nebula with the E VN and e-MERLIN com- 
bined would be a highly desirable option for a followup of any 
future strong flaring event in this object. 
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